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A GENERAL RELAXATION THEORY OF SIMPLE LIQUIDS* 


by Mati Merilo and E. J. Morgan 
CHAPTER I 
INTRODUCTION 

The liquid state is still the least understood of 
the three states of matter. The behavior of gases at 
low pressures can be well explained in terms of kinetic 
theory (Ref. 1). This is possible because the molecules 
of a dilute gas are so widely separated that they are 
unaffected by molecular forces, except, for very brief 
encounters or "collisions". 

Solids are affected strongly by intermolecular 
forces. However, because they possess an ordered struc- 
ture with the only motion of the molecules occurring as 
thermal vibrations, the mathematics required to analyze 
them is still tractable. 

Liquids,' due to their close spacing, are character- 
ized by strong interactions between molecules, yet they 
do not possess long range order, hence neither of the 
basic simplifications applicable to solids or gases can 
be used over the entire liquid range. 

Above the critical pressure, it becomes impossible 
to distinguish between the liquid and gaseous states . 
This has led to the approach of viewing liquids as dense 
gases, and much formal work on statistical mechanics and 

*Work sYppoTted by - the National Aeronautics and Space Administration 
under grant NGL 36-003-064 and the National Science Foundation under 
grant GK-02866. 



kinetic theory has been done from this viewpoint. 

(Refs. 1, 2, 3). 

Near the melting point, the molecules of a liquid 
exhibit close packing similar to the solid case. The 
difference between the two states concerns the molecular 
arrangement. In a crystalline solid the molecules 
vibrate about their centers, which are located in a 
regular geometrical pattern or space lattice. Therefore 
solids are said to possess long-range order. In a 
liquid, there is still a tendency for the molecules to 
exhibit ordering, however, this order does not persist 
beyond a few molecular diameters. This short range 
order has led to attempts to approach liquid theories 
from the solid side. Notable among these are the cell 
(Ref. 3) and hole (Ref. 4) theories. 

The dynamical behavior of fluids has been described 
with a good degree of success, for over a century, by 
the Navier-S t okes equations. However, when flow times 
become comparable to the times required for molecular 
processes, the equations are hopelessly inadequate. 

This problem was first encountered in measurements of 
absorption and dispersion of sound waves in gases. Good 
agreement between experiment and theory was finally 
provided by relaxation theory (Ref. 5) which was first 
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proposed by Herzfeld and Rice (Ref. 6). 

The basic premise of relaxation theory is that a 
finite time is required for excess energy in any one 
mode to redistribute itself among the other modes. Thus 
relaxation theory splits up the overall energy equation 
into separate equations for each mode and couples these 
equations by relaxation terms. Applying this concept to 
sound waves propagating through a gas, we see that energy 
is first imparted to the gas molecules as translational 
kinetic energy. After a number of collisions, this 
excess energy in the translational mode is redistributed 
among the rotational and vibrational modes. This process, 
known as thermal relaxation, leads to absorption in 
excess of that predicted by classical theory which con- 
siders only viscosity and conduction. 

Relaxation like behavior has also been observed in 
liquids, and in many instances the same relaxation 
mechanism as was found in gases has been found to apply 
(Refs. 7, 8). These liquids are sometimes known as 

Kneser liquids. A further type of thermal relaxation in 
liquids is caused by the perturbation of molecular equi- 
“ librium between different species , brought about by 
temperature variations. 
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In addition many liquids , primarily the associated 
liquids, display a relaxation caused by volume changes. 
The prime example of this is water, which shows excess 
absorption at 4°C. This absorption cannot be thermal in 
origin since at this temperature c^ = c^, and therefore 
adiabatic compression and rarefaction waves are also 
isothermal. For this reason the excess absorption in 
water has been attributed to the presence of both loose 
and close packed structural forms. The equilibrium 
amount of each structure is perturbed by changes in 
volume, and thus this process is referred to as a 
structural relaxation (Refs. 8, 9). 

Relaxation equations for liquids have also been 
written from purely phenomenological viewpoints to 
explain the viscoelastic behavior of substances. These 
are rheological equations of state which have been used 
to describe non-Newtonian flow, normal stress effects, 
and the propagation of ultrasonic shear waves. These 
theories are based on the suggestion by Maxwell (Ref. 10) 
that the behavior of a fluid can be described in terms of 
a relaxing solid. Oldroyd (Refs. 11, 12) has written 
general rules for formulating rheological equations of 
state, and has found that the simplest equations which 
are linear in stresses and include terms of second degree 
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in stresses and velocity gradients taken together, 
possess eight arbitrary constants. Though equations 
obtained from strictly phenomenological grounds can 
describe many observed phenomena, their main weakness 
is that there is no connection between the molecular 
properties of the fluid and the arbitrary constants 
required to describe the continuum behavior. This is 
in contrast to the relaxation theory of gases where a 
knowledge of molecular properties is used to obtain the 
relaxation equations. 

It is clear that any theory which hopes to describe 
the dynamical behavior of liquids must in some form / 
include relaxation effects. 

The purpose of this research is to relate the 
relaxational behavior of liquids to their molecular 
properties. The approach adopted is to evaluate the 
thermodynamic properties of the system by statistical 
mechanics, and to account for dissipative processes by 
relaxation from non equilibrium states towards equilib- 
rium. This approach is valid as long as the deviation 
of the radial distribution function from its equilibrium 
value is small (Ref. 13) , and this will be assumed to be 
the case in the following analysis. 
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CHAPTER II 


THEORETICAL ANALYSIS 

In the analysis that follows, two ideas are essen- 
tial and will receive further attention. One is that a 
liquid can, to a certain extent, support a recoverable 
shear strain, and thus a "molecular strain" tensor can 
and will be defined. This is possible because a finite 
time is required for a liquid to come to equilibrium 
following some disturbance, even though this time may be 
extremely short. 

The other important concept is that the random 
velocity, and thus the translational kinetic energy, of 
molecules can be a function of direction (Ref. 14). With 
a slight modification of the definition of temperature, 
this leads to the possibility of defining a "transla- 
tional temperature tensor". 

Both of these concepts will be used to write the 
stress tensor as a function of temperature and molecular 
spacing. Further, since both the strain and temperature 
will be related to the strain rate through relaxation 
equations, the stress willshow relaxation effects, and 
be a function of the strain rate. 
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2.1 Concept cf Strain in a Liquid, and the Molecular 

Strain Tensor 

A liquid, which is in static equilibrium with its 
surroundings has a uniform distribution of molecules 
throughout its volume. Thus we say that the radial 
distribution function, which specifies the number of 
molecules to be found at any distance from an arbitrary 
reference molecule, is spherically symmetric. If, as 
the result of some disturbance, the distribution function 
deviates from spherical symmetry, it will take a finite 
time for molecular collisions to bring the distribution 
back to its equilibrium value. Another way of viewing 
this is to note that as a result of some disturbance the 
average spacing between molecules can vary with direc- 
tion . 

It is possible to draw some analogies between a 
liquid with a distorted radial distribution function and 
a strained solid. In discussing solids. Love (Ref. 15) 
states that "the state of strain is entirely determined 
when we know the lengths in strained and unstrained 
states of corresponding lines". Thus for a solid, in 
principal axes of a cartesian system, a strain tensor can 
be defined in terms of extensions of material line 
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elements, or fibers, along the principal directions. If 

s is the initial unstrained length of a line in the a 
OCX 

direction, and s q is the strained length, the elongation 
of the fibre or line is defined as 


Y 


a 



( 2 . 1 ) 


There is no unique measure of strain (Ref. 16) , and in 

fact any function of Y a can be used provided that it 

vanishes for y =1, reduces to the infinitesimal strain 

tensor for |y — 1 I < < 1 , and is dimensionless. As a 
1 'a 1 

result, many measures of strain appear in the literature. 
Some of the more common ones are (Ref. 16): 


Cauchy measure 


e = Y - 1 
aa 1 a 


Hen cky measure 
H 


e = Jin y 
aa a 


Green measure 



e aa 2 * . . 

plus many others associated with the names of Swainger, 
Alamans i et c . . . . 
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For a liquid, if the length of the material line, 
or fiber, is replaced by the average distance between 
molecules in a particular direction, it is possible to 
describe the distortion of the radial distribution 
function in terms of the changes in molecular spacing 
along principal axes. We can then define a "molecular 
strain" tensor in principal axes, a, as 



( 2 . 2 ) 


where £^ is the average distance between molecules in 
the a direction and £ is the average molecular spacing 
at the local density. Since this tensor is expressed in 
principal axes of the deformed material, it can be 
related to the Cauchy deformation tensor (Ref. 17), ^c, 
by 

e + 1 = c ^ 


where £ is the molecular strain tensor, and ^ is the unit 
tensor. The tensor ^ ^ is convenient because the prin- 
cipal values are exactly the principal elongations, how- 
ever, in axes other than principal axes, the terms of the 
tensor are considerably more difficult to evaluate in 
terms of the deformation gradients. In fact, the compo- 
nents are in general complicated infinite series in 
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displacement gradients (Ref. 17).- The change in volume 
during a deformation is expressed by 


d V 


-h 

det (c) dV 


Since the molecular strain refers only to distortion, 
and not to dilatation the condition 

det (6* + ej) - 1 (2.3) 

must be satisfied. 


2 . 2 The Stress Tensor 

From the kinetic theory of dilute gases, the pres- 
sure tensor can be expressed as 

p.. = n ( m w . w , } 

ij ' i y 

where the subscripts refer to the axes of a cartesian 
coordinate system, n is the number density, m is the mass 
of the molecule, and w is the random, or peculiar, 
velocity. The usual kinetic theory definition of 
temperature is just 

(2 - 4 > 

where k is the Boltzmann constant. 
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However, this definition can be extended such that a 
translational temperature tensor can be defined as 




w 


>* 


(2.5) 


This leads to an equation of state for an ideal gas 
(Ref. 14) 


P 


ij 


p R T 


ij 


( 2 . 6 ) 


where p is the density and R is the gas constant. The 
traditional temperature then is just one third of the 
first invariant of the temperature tensor. 

Similar equations of state can be written for a 
liquid, incorporating the effects of the attractive 
forces and the finite size of the molecules. 

Since directionally dependent molecular spacings 
and temperatures must, if not acted upon by external 
forces, eventually equalize, we see that we are in fact 
dealing with a non-equilibrium situation. The usual 
approach in studying non-equilibrium statistical 
mechanics is to investigate the change of the distribu- 
tion function towards canonical form as a solution of the 
time dependent Liouville equation. The main problem, so 
far, has been to calculate the distortion of the equilib- 
rium radial distribution function for specific non- 
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uniformities. Due to the complexity of methods based on 
the Liouville equation, there is not yet available a 
complete theory of simple non-uniform liquids that leads 
to results which can unambiguously be compared with 
experiment. With the hope of achieving mathematical 
tractability this approach is abandoned in favor of a 
simpler, more intuitive procedure. To describe the 
directionally dependent temperatures, Morgan and Kern 
(Ref. 14) hypothesized that the non-equilibrium distribu- 
tion function could be specified by a quasi-Maxwellian 
distribution which has the property that it is an 
equilibrium distribution in each direction taken sep- 
arately. In this way it is possible to talk about 
thermodynamic quantities as functions of direction. From 
this viewpoint, the irreversibility of the flow arises 
from the coupling of the different directional modes by 
relaxation terms, and allowing the system to return 
towards the isotropic equilibrium distribution. 

The results obtained by Morgan and Kern, for the 
directional dependence of temperature, are used here. 
Further, we postulate a quasi-equilibrium molecular 
spacing in the principal axes of deformation, an assump- 
tion which allows us to use the powerful techniques of 
equilibrium statistical mechanics. There are two 
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distinct but equivalent procedures which may be applied. 
One calculates the thermodynamic properties through the 
partition function, and the other through the radial 
distribution function. The former approach will be used 
here . 


There are a number of ensembles which are used for 
thermodynamic calculations, depending on what the inde- 
pendent variables are, and the types of contact the 
systems make with the surroundings. It can be shown 
(Ref. 18), however, that away from critical regions, 
where fluctuations in thermodynamic properties are small 
the ensemble to be used can be chosen for mathematical 
convenience. Generally, the simplest and therefore most 
commonly used one is the canonical ensemble. 

The canonical partition function is 


Z n (V,T) 


I exp 
i 



(2.7) 


where N is the number of molecules, V is the volume, and 

\ 

is the energy of each system of the ensemble. The 
summation is taken over all systems of the ensemble. 

Once the partition function is known, it is possible to 
obtain all the thermodynamic properties from it-. r 'In 
particular, the Helmholtz free energy is given by 

F = - kUn Z . (2.8) 
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From this, the thermodynamic pressure can be expressed 
as 


P = - 



(2.9) 


The term pressure will only be used for the case of an 
unstrained liquid. When the liquid is strained, since 
the molecular strain is recoverable, elasticity theory 
can be used. Truesdell (Ref. 19) has pointed out that 
there really is no satisfactory thermodynamical treatment 
of the foundation of elasticity theory, therefore it is 
necessary to postulate that a stress tensor obtained 
from an "elastic potential" or "strain energy" is in fact 
equivalent to the stress tensor obtained from the prin- 
ciples of mechanics. When this is done, the elastic 
potential for an isothermal deformation is given by the 
Helmholtz free energy (Ref. 19), and the stress tensor, 


Sj , can be expressed as 


S* = I x 1 3F 


J V ’ K 3xJ 


( 2 . 10 ) 


,K 


Here x is the displacement or deformation gradient 
» R 


defined by 


i _ 3x 


,K ax K 


where X and x are coordinates in the undeformed and 
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defonned material respectively. In principal axes we can 
relate the molecular spacing to the displacement gradient 
with the expression 

x ct = 

*“ ax a a a 

o 

where Jl is the equilibrium spacing in some arbitrary 
reference state. Thus it is possible to relate the 
stress tensor, through the Helmholtz free energy and the 
partition function, to the molecular properties of the 
liquid . 

2 . 3 The Momentum Equation 


The momentum equation is used in its usual form 

.i 



Du 1 
P Dt 

where u 1 is 

the ve 

nent of any 

body f 

refers to c 

ovarian 

substantial 

deriva 


Da' 1 ; ' 
Dt' 0 ' 


Pf‘ 


( 2 . 11 ) 


3a' i ; ’ k .i. . 
= ^..j. + u a> 


( 2 . 12 ) 
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2 . 4 The Molecular Strain Relaxation Equations 


Equations which express relationships between the 
stress tensor and the kinematic variables, such as 
velocity and acceleration, at any point of the body are 

called rheological equations of state. Since the stress 

\ 

tensor can be expressed in terms of the molecular strain 
tensor and the temperature tensor, the relationships 
connecting these tensors to the strain rate tensor can 
properly be called rheological equations of state. These 
equations describe the properties of an arbitrary material 
element moving as part of a continuum and as such must be 
independent of the frame of reference. 

An arbitrary fluid element in motion is continuously 
translating, rotating and deforming, however, the stress 
in the element should only be affected by the strain. If 
this were not the case, the stress distribution would be 
a function of the rigid body motion even if inertia and 
body forces were neglected. To avoid the effects of 
rigid body motion, the physical behavior should be 
expressed in a coordinate system which is rotating and 
translating with the element. The only motion which 
would be seen in this system is a pure deformation, and 
thus strain and relaxation effects could be more easily 
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expressed. It is not very convenient to solve problems 
in this system, however, thus the results should be 
transformed into standard Eulerian coordinates. To sum- 
marize, the procedure to follow is to state a physical 
law, or hypothesis, in a coordinate system following a 
fluid element, and then to transform this to a fixed 
coordinate system for solving problems. 

This transformation is easily accomplished by 
replacing the time derivatives appearing in the moving 
system by the Jaumann derivative defined as (Refs. 12, 
20 ) 


Da ' 
D t ' 


• J • 


i. . 







(2.13) 


where \ denotes that the summation over m must be per- 

_ » 

formed for each covariant index and l denotes the same 
for each cont ravariant index. The vorticity is defined 
in the usual manner as 


0) 




U j,i ) 


(2.14) 


The Jaumann derivative can also be written in terms of 
the more familiar material derivative-. It is - 


Da ' * 
Dt ' • 


Da . ” ' r m 

= 5F- 1 ** “ l w i a . 


• . j 


• j . 

m . . 


r i . . m , 
-> a . 

L m . i . , 


(2.15) 
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where the extra terms account^ for the rotation of the 
fluid element. It should be noted that for a scalar 
quantity the last two terms of equation (2.15) vanish, 
and the Jaumann derivative reduces to the common material 
derivative. 

The strain of a body can be described by the rela- 
tive motion of two points when the rigid body motions 
have been removed. This motion can be considered to 
consist of three strains along the principal axes of the 
strain-rate tensor. Thus it can be shown (Ref. 21) that 
if 6s is the distance between the points, its relation- 
ship with the strain-rate tensor is just 


-r— a = d (no summation) (2.16) 

6s D t aa 

a 

where d is the aa component of the strain-rate tensor 
aa . . 

defined as 



7 <u i.J + u J,i ) 


(2.17) 


This equation, as it stands, shows that a medium with a 
constant strain rate imposed on it can deform, or strain, 
indefinitely. However, the molecular strain, except for 
the case of pure dilatation can not. This is because the 
molecular strain depends only upon the mean configuration 
of molecules about a given reference molecule, irrespec- 
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tive of the identity of the surrounding molecules. Thus 
when two molecules are separating from each other because 
of a strain-rate, and another molecule moves between 
them, the strain is effectively relieved since the new 
neighbor determines the radial distribution. Upon cessa- 
tion of the strain-rate, the radial distribution function 
will after a finite time revert back to its spherical 
symmetry, and thus the molecular strain will relax to 
zero . 

If in equation (2.16) the length Ss^ is replaced by 
the average distance between molecules in the a direc- 
tion, £ , extra terms must be added to account for the 
relaxation behavior. We shall postulate that the relaxa- 
tion can be expressed by terms which are linear in the 
deviation from equilibrium. We further make the assump- 
tion that the molecular spacing in one direction can be 
a function not only of the strain rates in the same 
direction, but also of those normal to it. Then, in a 
cartesian system, in principal axes of the strain, we 
can write 


1_ 

t 


a 


Dl , 

— a + I 

D t < 


fi -a 

a x 


a -a 

_« I 


i -l 
a z 


V 


= (A-B)d + 

aa 


B ( d + d 
xx yy 


+ 


d ) 
z z 


(2.18) 
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where K is the relaxation time for the equilibration of 

molecular strain and A and B are coefficients for the 

effects on molecular spacing caused by direct and normal 

strain rates respectively. One would expect that for 

simple liquids cross strain effects should be small, thus 

A should be much larger than B. We see that when 

l = l = l , the relaxation term is zero. However, when 
x y z 

there is a difference in spacing in different directions, 
relaxation will occur to decrease this inequality. In 
terms of the molecular strain tensor, equation (2.18) can 
be expressed as 


.a D£_ 
’a Dt 


aa 


aa k 

+ U-.)d“ + ,dJ| + 3 l§| 


(2.19) 


where is the inverse of (e^ +<$^). It is evaluated 

J i i 

from the expression 


c 1 - 

J e 


( 2 . 20 ) 


where is the cofactor of the element ej + 6*? in the 
j IX 

determinant (e^ + 6^) = £. In an arbitrary coordinate 
system equation (2.19) becomes 


1 M oe kj + ss K1 1 

2 r k Dt Dt ) 


ki i 3 e 
+ — 


ij - g ij e k ii 

8 __k _ &_L DP = 


K 


3p Dt 


= (A-B)d 1J + Bg ±j d k , 


( 2 . 21 ) 
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ana if we lower j , set i = j and contract we get 

- ^ = (A + 2B) d£ (2.22) 

£ Dt p Dt k 

But from equation (2.3) we have £ = 1, therefore 
equation (2.22) becomes 

p Dt + (A + 2B) d k = ° * (2.23) 

and for this to satisfy continuity we must impose the 
condition 

A + 2B = 1 . (2.24) 

Since we expect that A >> B, A should be almost unity 
and B should be close to zero. 

2 . 5 The Energy Relaxation Equations 

If in a stationary gas the random translational 
kinetic energies of the molecules are unequal in dif- 
ferent directions, molecular collisions would eventually 
redistribute the energy until the inequality vanished. 
Differences in the temperatures would thus relax to zero. 

Morgan and Kern (Ref. 14) have shown from a kinetic 
theory analysis that in principal axes of the strain 
rate, separate energy relaxation equations can be written, 
for a gas, in each of the principal directions. In a 
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cartesian coordinate system where a represents any prin- 
cipal direction, their expression with their simplified 
conduction terms is 


pc ( DT 


aa 


(T - T ) + (T - T ) + (T 
aa xx aa yy aa 


Dt 



( 2 . 25 ) 


where t is the relaxation time for energy exchange 
between translational modes, c is the specific heat at 
constant volume, and k is the thermal conductivity. It 
should be noted that if the three equations in the three 
principal directions are added, the usual energy equa- 
tion for an ideal gas is recovered. The form of the 
energy equation for this case is contingent on the fact 
that the internal energy for an ideal gas is a function 
of temperature alone. If, as is the case for a liquid, 
the internal energy depends also upon the density and 
the strain, the form of the energy equation is altered 
by the addition of several new terms. 

The usual form of the energy equation is 


A 

p Ilf ' + + si3 "j,i (2 - 26) 

A 

where E is the internal energy per unit mass. 
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the left-hand side of equation (2.26) 


For E = E(p,T,e ) , 
can be expanded as 
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Dt 
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L 3 pj 

T ' £ iJ 
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T Dt 

If we make the assumption that for small strains 

f 3El 


(2.27) 
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ij 


f 3E' 
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and 


(2.28) 
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Utj 

k 

p. 


(2.29) 


and if we define II ^ as 


il ^ P(«, + 


3E 


3e 


P,T 


(2.30) 


then, if we also make use of the continuity equation the 
left-hand side of equation (2.26) can be rewritten as 


D§ DT , 

P Dt PC v Dt + 


'--fir) 


3T 

‘ ' P j 


d k + " lk 57“ 


C2. 31) 


Equation (2.30) is seen to define a stress tensor whose 


elastic potential is just the internal energy. 
If we further define the tensor by 
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s lj = -Pg 1 - 3 + 


C2.32) 


where P is the thermodynamic pressure for the unstrained 
material at the local mean internal energy and density, 
the energy equation becomes 


T l + 0% 5T - + < £ e Imi , 4 ).. + oi3d 


ji 


(2.33) 


where k is the thermal conductivity and 


M _ £ 
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f 3E 


De: 3 
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/ ik J 

P,T 


^ jk ; 

P,T 


(2.34) 


is the rate of change of strain energy. Except for the 
strain energy term, equation (2.33) is a common form of 
the energy equation (Ref. 21). We now assume that it is 
possible to write energy relaxation equations for a liquid 
which are analogous to equation (2.25) derived from 
kinetic theory by Morgan and Kern. We propose the 
equation 



+ 


3T 1:5 



T 





+ 


+ m + + ° Jk *i) < 2 - 35 > 


which reduces to equation (2.25) for an ideal gas, and 
can be transformed into equation (2.33) by lowering j, 
setting i = j, and contracting. 
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It is sometimes convenient, particularly for 
isothermal flows, or when conduction can be neglected, 


to use a non-dimensional deviatoric translational 
temperature tensor, 0^, defined by 


ij 


T ij ~ 8 lj T e_ 


(2.36) 


where 


T = 


I T k 

3 k 


(2.37) 


With this definition we see immediately that 


0 k = 0 (2.38) 

i 

This tensor is useful in that it expresses deviations of 
temperatures in different directions, rather than the 
actual temperatures. As such it is analogous to the 
strain tensor which expresses deviations from the 
unstrained conditions. Using this definition the energy 
relaxation equations can be rewritten as 


pc T 
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(2.39) 
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It should be noted that these equations do not consider 
vibrational and rotational energy modes, and should 
transfer between these and the translational modes be 
important, relaxation equations connecting these modes 
should be used . 

We have in this chapter expressed twenty-one 
equations which, in principle at least, can be solved 
for the twenty-one unknowns describing a general flow 
problem. The equations consist of six molecular strain 
relaxation equations, six energy relaxation equations, 
six stress equations and three momentum equations. The 
unknowns are six components each of the molecular strain 
the temperature, and the stress tensors and three 
components of velocity. Before a solution to these 
equations can be obtained, however, it is necessary to 
evaluate the partition function, and to do this a model 
for the liquid must be chosen. This will be done in the 
next chapter where the cell theory approach will be used 
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CHAPTER III 


THE CELL THEORY APPROXIMATION 

To obtain solutions to the equations derived in 
Chapter 2, it is necessary to evaluate the partition 
function for a large system of interacting molecules, a 
task which at present is quite unfeasible. This diffi- 
culty is what led Eyring and Hirschfelder (Ref. 22), and 
Lennard- Jones and Devonshire (Ref. 23) to attempt to 
simplify the evaluation of the partition function by 
introducing cell theory. Though the theory was orig- 
inally proposed as a result of an intuitive approach, 
Kirkwood (Ref, 24) has placed the theory on a sound 
statistical mechanical basis by obtaining the partition 
function for the cell model by successive approximations 
of the exact partition function. 

3 . 1 The Basis of Cell Theory 

The physical arguments leading to the adoption of 
cell theory emphasize the similarity of the liquid struc- 
ture to that of a crystalline solid. In a crystal," the 
molecules are located in a regular geometrical pattern, 
and therefore the crystalline solid is said to exhibit 
long range order. In a liquid, for high densities and 
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for temperatures far below the critical point, X-ray 
diffraction studies have shown that short range order, 
extending over a few molecular diameters, still exists. 
Since the mean molecular spacing is not much greater than 
the diameter of the molecule, it can be assumed that a 
molecule is confined to a "cell" or cage consisting of 
its immediate neighbors. This is a reasonable assumption 
for thermodynamic calculations provided that the mole- 
cules occupy their cells for time periods which are large 
compared to the time between collisions. A molecule 
which is restricted to move in its cell, can thus be 
regarded as an independent thermodynamic system. If it 
is assumed that the system is homogeneous in the sense 
that each molecule in a cell is representative of all 
other molecule-cell combinations, then the partition 
function for a system of N molecules can be obtained by 
taking the Nth power of the single molecule partition 
function. In terms of the Helmholtz free energy, this 
just means that the total free energy is N times the 
free energy of one molecule in one cell. 

In spite of the number of approximations, the cell 
theories have had outstanding success in obtaining 
qualitative and, to a lesser extent, quantitative agree- 
ment with experimental thermodynamic properties. It is 
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possible to achieve varying degrees of complexity with 
cell theories, depending on whether the number of mole- 
cules is put equal to the number of cells, whether 
single or multiple occupancy of the cells is permitted, 
and depending on the form of the intermolecular potential 
which is chosen. 

3 . 2 Evaluation of the Partition Function and the 
Thermodynamic Properties 

To facilitate the mathematics, the simplest form of 
cell theory will be used here. We assume that there are 
N molecules occupying N cells, with only single occupancy 
permitted. Attention will be focused on only one mole- 
cule, which is allowed to move, while the remaining mole- 
cules are assumed to be fixed in their equilibrium posi- 
tions in a quasi-crystalline lattice. The free molecule 
or the "wanderer" is restricted to move within its cell, 
where the potential is assumed to be uniform. Thus the 
potential energy can be expressed as 

<j> = — (within the cell) 

<j> = °° (outside the cell) 

where $ is the total lattice potential energy of the 
Li 


L,. (3.1)-.. 
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molecules in their equilibrium positions. It is gener- 
ally expressed as a function of the specific volume. 

The canonical partition function, defined by 
equation (2.7), can be written as a product of terms 
including the translational kinetic energy and the 
potential or configurational energy. Thus 


where 


Z = Z 


kinetic 


Z 


kinetic 


[ 2’H’ mkT 1 


3N 

2 


) 


(3.2) 


(3.3) 


h is Planck's constant. The conf igurational partition 
function for N particles is 


Q n (V,T) 



. . . . r^) 
kT 


dr , . . . . dr M 
—1 — N 


(3.4) 


where r. is the radius vector of the i ' th molecule, 

—l 

$(r^ . . . . _r N ) is the potential energy of the system in 
configuration (r^ * • • •— ^) and the integrations are carried 
out over the volume V. Utilizing the cell theory assump- 
tions, the configurational partition function for one 
molecule in its cell is 


Q 1 




dr 


(3.5) 
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where r is the position vector of the molecule in its 


cell. Since 4> is assumed to be uniform 


Q 1 = v fc exp 


f- M 

l kTj 


(3.6) 


where v, is the free volume available to the center of 
f c 

the wanderer in its cell. The configurational free 
energy for one cell is then 

F, = <b - kT An v r C3.7) 

1 f c 

and for N particles in N cells 

F„ = $ T - NkT An v _ (3.8) 

N L tC 

If it is assumed that the immediate neighbors 
making up the cell are fixed in their lattice positions, 
then for a face centered cubic lattice the volume per 
molecule is 

v - A 3 //2 (3.9) 

where A is the distance between nearest neighbors. The 
free volume available for the center of the moving 
molecule has a very complicated shape, however, the error 
introduced is not large if the volume of the largest 
sphere which fits inside the actual free volume is taken 
as the free volume (Ref. 25). 
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Then 



y it (£,-d) 3 


(.3.10) 


where d is the hard sphere diameter. When the lattice 
is strained, the relative distances between the molecules 
and the shape of the free volume are altered. Choosing 
cartesian axes x, y, and z, coincident with the principal 
axes of strain for the lattice, the sphere representing 
the free volume is distorted into an ellipsoid of volume 


v c = — 7T X X X 

fc 3 xx yy zz 


(3.11) 


where 


X = (1 + e H -d . 
aa ' aa e 


(3.12) 


Thus X represents the size of the gap between the 

(X \X 

molecules in the a direction. 

To evaluate the partition function, it is necessary 
to choose an expression for the interaction potential 
between molecules. There are many forms of the inter- 
molecular potential currently in use, giving good descrip 
tions of gas viscosities, second and third virial coeffi- 
cients, thermal conductivities and other properties. 
Unfortunately, since the functions used are only empir- 
ical representations of the true interactions, the 
parameters required to describe them depend on the 
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property being measured. As an example, the parameters 
of the Leonard- Jones (6-12) potential required to describe 
gas viscosities are different from those required for the 
second virial coefficients. At present, there does not 
appear to be any potential function which allows consist- 
ent information about intermo lecular properties to be 
derived from various physical properties. In view of 
these uncertainties, a simple form of the potential will 
be chosen, which does, nevertheless, represent to some 
extent the actual behavior. 

It has been shown (Ref. 3) that a good approximation 
to the lattice potential energy is given by a function of 
the form 


<j> = - (3.13) 

n 
v 

For a liquid which is being strained, . plausible expres- 
sion which reduces to the above form for no strain is 


where 



• 1 1 1 ' 

^ + 7 ^ + 7 ^ 

l x y z 


2 n/2 


a 


(3.14) 


(3.15) 


and 1^, Z ^ represent the molecular spacing in the 

principal x, y, and z directions respectively. It is now 
possible to evaluate the desired thermodynamic properties. 
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The stress tensor can be obtained in principal axes as 
follows , 


aa 


"a [~3F 
V 95, 


C3.16) 
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Substituting for F from equation (3.8) we get 
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Z kT 

'9 

In v £ 
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(3.17) 


where the effect of temperature variation with direction 
has been included. When the derivatives are evaluated, 
the stress becomes simply 


S = - pRT 
aa aa 


1 + 


aa 


na 


n+1 j. ' - 3n 
v (1 + e ) 
aa 


(3.18) 


For the case of zero strain, and n = 1, we get 


- S = P = pRT 1 + f 
aa X 


(3.19) 


which is one form of the Eyring equation of state 
(Ref. 4). The identical expression for pressure can be 
obtained by assuming no strain and using the relation 



(3.20) 


A closer look at equation (3.19) indicates that in the 
low density limit, the equation of state for an ideal gas 
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is recovered. The reason for this behavior is that for a 


system of non-interacting particles, the introduction of 
partitions affects the free energy and the entropy, but 
not the pressure or the internal energy (Ref. 3). Thus 
equation (3.19) is correct in the high and low density 
limits, though it can be expected to be in error in 
between . 

The stress and the pressure may be expressed as the 
sum of kinetic and potential terms. The kinetic contribu 
tion is a consequence of the motion of the molecules in 
their cells, while the' potential term results from the 
potential energy of the molecules located at their respec 
tive cell centers. Thus the kinetic pressure is 


P 


k 


pRT 1 



(3.21) 


and the potential pressure is 

p p'-Si ‘ 3 - 22) 

v 

For small strains, the potential term of the stress 
tensor can be linearized so that equation (3.18) becomes 


S = - pRT 1 + 


ota 


aa 


aa 


+ P (l-3ne ) . 
p aa 


(3.23) 


We now define a non-dimensional gap tensor,, , as 
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(3.24) 



+ P p (g ij -3ne ij ) (3.26) 

It is often more convenient to use the deviatoric 
stress tensor defined in equation (2.32) and to write it 
in terms of the deviatoric translational temperature 
tensor defined in equation (2.36). Equation (3.26) can 
then be rewritten as 

-1 

°ij ‘ - pRT e( 0 iJ + ilKk + 8 ik ] V J + 

-l 

+ t 0 jk + V 1 ^ ' 2 e U }) ' 3nP P e ij • (3 ' 27) 
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To find an expression for 11^ , which was defined by 
equation (2.30), we must first determine the internal 
energy, E. This is most easily accomplished by using 
the Gibbs-Helmholtz equation 


E = F - T 



(3.28) 


From the partition function which was obtained 
earlier, the free energy is 


F = - iSM m 


2irmkT 


- NkT £n (v, ) + 

X C L 


(3.29) 


and the internal energy becomes simply 


E = f NkT + $ L 


= c T + $ T 
v L 


(3.30) 


Using the definition for 11^ , we have in principal 
cartesian axes that 


n 


act 




T,P 


(3.31) 


which turns out to be simply the potential part of the 
stress tensor 




3ne 



(3.32) 
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and for no strain 



3.3 The Relaxation Times 


(3.33) 


In the preceding chapter, two relaxation times were 
postulated, one for equilibration of translational energy 
and the other for equilibration of molecular spacing, but 
no indication was given of how they may be obtained. The 
time required for equipar tition of translational energy 
for molecules of equal mass should be of the order of the 
collision time. Morgan and Kern (Ref. 14) have shown 
that in a gas the translational relaxation time is just 


T 



(3.34) 


thus if the kinetic theory result for the viscosity of a 
hard sphere gas is used, it is found that 


T = 3.813 T (3.35) 

c 

where T is the collision time. It will be assumed that 
c 

the same relationship between the relaxation and collision 
times will apply for a hard sphere liquid. 

The collision time for a liquid is a very difficult 
quantity to determine. Herzfeld and Litovitz (Ref. 5) 
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have shown that collision times computed under various 
assumptions for some common liquids disagree by factors 
of up to thirty. In view of this uncertainty, the col- 
lision time will be calculated using a convenient though 
relatively simple approach. The average collision fre- 
quency for molecules in a liquid can roughly be expressed 
as 


C x C C z 

f c = 2X + 2X + 2X 
x y z 


(3.36) 


where c^ is the component of the average random speed of 
the molecules in the a direction. By the component of 
the random speed we mean the component of the absolute 
value of the random velocity. Since the velocity distri- 
bution in a liquid, as in a gas is Maxwellian (Ref. 33), 
the average speed is 



(3.37) 


At equilibrium, the average speed is independent of 
direction, thus the component of the speed in any direc- 
tion is obtained by averaging the component for the 
velocity vector varying over all angles, thus 

c n = ~ W (3.38) 

0 7T 2 

The spacing is also uniform at equilibrium, thus the 
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collision time can be expressed as 


T 

C 


1 it 2 \ /irn 

f = 12 J 2kT 

c * 


(3.39) 


When the liquid is strained, the spacing between 
the molecules in different directions is altered and, as 
may be seen from equation (3.36), the collision frequency 
is changed. 

The ratio of the equilibrium to the strained 
collision frequency is 


f 3c/ A 

ce 



(3.40) 


It will be shown later that 
therefore 0.. << e . . < V . 

ij ij ij 

vary as the square root of 
of strain on the collision 
effects of the temperature 
assumed that the molecular 
value even in the strained 


t 

in most liquids T << K , and 
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the temperature, the effect 
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In any axes 
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where the invariants of the gap tensor are 

IX = — v ^ - v* vM 

V 2[ V i V j V j ij 


and 


III = det 
V 


ft) 


(3.43) 


(3.44) 


The structural relaxation process is considerably 
more complicated than the thermal one. When a material 
is strained, the separation between molecules is altered 
The size of these separations fluctuates continuously 
such that at some time a vacancy large enough to accom- 
modate a molecule is formed. When a molecule jumps into 
this available space the strain is relieved. Therefore 
it is possible to identify the relaxation time for mole- 
cular strain with the time required for a molecular jump 
Two conditions are required for a jump to occur. 
First, there must be a sufficiently large volume avail- 
able, and second, the molecule must have enough energy t 
leave its position and surmount the potential barrier 
formed by the neighboring molecules. 

Chung (Ref. 26), using statistical mechanical argu- 
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Brummer (Ref. 28), and Gubbins and Tham (Ref. 29) 
have shown that the energy required for a jump is a strong 
function of the density of the fluid. This is to be 
expected since a molecule must squeeze through the 
extremely steep repulsive potential of its neighbors. In 
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a strained liquid then, where a molecule sees a density 
varying with direction, the total jump probability is 
obtained by integrating over all directions. Thus 


p j 



+ E v('K'» 

kT 



E v(<t>,iJO 

kT 


v f C 4> > ^ ) J 


sin<j)d<J>diJ; 


C3.48) 


where (j) and iJj are the angles measured in spherical 
coordinates . 


3 . 4 Limitations of the Simple Cell Theory 

Although cell theory has been shown to give a useful 
description of the liquid state, there are nevertheless 
several limitations introduced by the many assumptions 
which are made. 

One of the main problems is that cell theory is 
really a description of the solid state, introducing a 
degree of long range order which does not exist in a 
liquid. The introduction of cells implies that mole- 
cules in one cell do not affect molecules in neighboring 
cells and thus correlation effects are neglected. Both 
of the above assumptions lead to predictions of entropy 
which are too low, since order has artificially been 
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introduced to the system. 

The supposition that the virtual lattice, defining 
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The range of applicability of cell theory is limited 
by the specification of the cell center as the preferred 
position for a molecule. Calculations for molecules with 
a Lennard-Jones (6-12) potential indicate that at reduced 
density 

; 

p* 5 nd 3 < 0.69 (3.49) 

where n is the molecular number density,' there is a 
potential maximum at the center of the cell CRef. 32). 
Since for this case the ground state no longer corres- 
ponds to a regular arrangement of molecules, it is 
difficult in this range to justify cell theory. If the 
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energy of the molecule is less than the energy of the 
central hump, the effect would be for molecules to 
cluster, and in a sense behave like larger more compli- 
cated molecules. 

A common problem for theories of the liquid state 
arises when the potential contributions to the pressure 
become of the same order as the thermal or kinetic terms. 
Then it becomes necessary to subtract large terms of 
equal magnitude, leading to a magnification of any errors 
Physically the effects of this delicate balance are shown 
dramatically by the phenomenon of critical opalescence, 
an occurrence which is brought about by large scale 
fluctuations in density. The density fluctuations are 
in turn a manifestation of the formation of clusters by 
attractive forces and their subsequent dispersion by mole 
cular collisions. The clusters are found to contain of 
the order of millions of molecules (Ref. 33). Clustering 
to a lesser extent must occur in liquids even away from 
the critical point, thus, in particular, near the satura- 
tion curve the cell theory cannot be expected to be 
accurate . 

’ The limit a t ion s - imposed by the ad option of ce 1-1 - 
theory must be kept in mind when interpreting the 
solutions obtained for various problems. 
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CHAPTER IV 


SOME SOLUTIONS AND COMPARISONS WITH EXPERIMENTS 

The problem with the present equations, as with most 
rheological equations of state, is that they are too com- 
plex to solve for any but the simplest of flows. The 
general approach in problems of rheology is to postulate 
a flow pattern, and then attempt to calculate the stress 
pattern necessary to maintain this flow. This procedure 
is used because the inverse problem of calculating the 
flow from the stress pattern, while difficult for 
Newtonian flow, becomes quite impracticable for more 
complicated flows. 


4 . 1 Steady Shearing Flow at Low Strain Rates 


Consider a flow such that in cartesian coordinates 
the velocity field is given by 


u 

x 


Dy 


u = u = 0 
y z 


( 4 . 1 ) 


where D is called the rate of shear. The coordinate 
system is shown in figure (1). 

Since the flow is described in cartesian coordinates. 
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it is no longer necessary to differentiate between the 
covariant and contravariant components of the tensors. 
Also, the metric tensor can just be replaced by the 
Kronecker delta. 

For this flow, the strain rate and vorticity tensors 
can be represented by 


and 



0 D/2 0 

D/2 0 0 

0 0 0 



0 D/2 0 

-D/2 0 0 

0 0 0 


(4.2) 


(4.3) 


respectively. The flow is further assumed to be steady 


3_ 

at 


o 


(4.4) 


isochoric (constant volume) 


d 


kk 


0 , 


(4.5) 


and isothermal. 

It is instructive to solve the equations at first 
for the case 


D 0 , 
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since then the terms of the strain and deviatoric temper- 
ature tensors are small, and the relaxation equations can 
be linearized. 

Equation (2.21), describing the relaxation of mole- 
cular strain, becomes 


^ £ i1 ^e. . ” ^il £ kk 

_ 3 + 1 . 1 - -A J B d kk . (4.6) 

since for low shear rates, and thus small strains 
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In terms of the individual components we can write for 
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for i = y, j = y 

(e - e ) + (e - e ) 
— yy zz 


d e + 
xy 
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(4.8) 


for i = z, j = 


(e - e ) + (e - e ) = 0 , 
zz xx zz yy 


(4.9) 


and for i= x, j = y 


? (e - e ) + = ( a-b ) ? 

2 xx yy' K 2 


(4.10) 
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The z axis remains a principal axis, therefore the 
off-diagonal terms with a z component are all zero. 

When the expression obtained from equation (4.9) for £ 

z z 

is substituted into equations (4.7) and (4.8), the 
latter become identical. This happens because the mol- 
ecular strains are not independent, but must satisfy 
equation (2.3). For small strains this gives 


£ 2 = e + e + £ . (4.11) 

xy xx yy zz 

Solving equations (4.7), (4.9), (4.10) and (4.11), the 

components of the molecular strain tensor can be 
expressed by the array 




(4.12) 


We see that the magnitude of the molecular strain is 
determined by the dimensionless product <D, a number 
which essentially expresses the ratio of the time 
required fox mo lex u Tar rearrangement to the~time "scale 
of the imposed fluid deformation. Using this number we 
are in a position to state that a requirement for a 
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shear rate to be considered- s-low Is that the condition 
KD << 1 must be fulfilled. For most simple liquids under 
normal conditions, K is extremely small, generally being 

_ g 

less than 10 seconds. 

Another important point to note is that the shear 
strain is of order <D, while the normal strains are of 
order (kD) 2 . This means that the expression obtained for 
volume changes in linear theory (Ref. 17), 


dV - dV 


(4.13) 


o 

cannot be used even for small strains, since e is of 

xy 

the same order as £. , . This is the reason for the £ 2 

kk xy 

term appearing in equation (4.11). 

The linearized energy relaxation equation can be 
written as 


Pc 
v_ 

3 





(4.14) 


where it has been assumed that the dissipation term is 
of order D 2 , and can thus be neglected. This assumption 
is verified a posteriori. In terms of the individual 
components, equation (4.14) can be expressed as. 
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for i = x, j = x 


XX 


IA 

3 9 xy 


( p \ 


= + 


pc T 
v 


tD e , 
xy 


(4.15) 


for i = y, j = y 


e + ~ 0 

yy 3 xy 


pc T 
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TD £ , 

xy 


(4.16) 


for i = z, j = z 


0 

z z 


0 , 


(4.17) 


and for i 


x, j = y 


0 + ^ (0 -0 ) = - 
xy 6 xx yy 


{ P 


pc T 
v 


tD , x 

-r- ( £ - £ ) - 

2 xx yy 


f 9P) 


P c , 


9T 


TD 

2 


(4.18) 

/ 


The equations for i ^ j and either i or j = z, become 
identically zero since z remains a principal axis. From 
equations (2.38) and (4.17) we find that 


0 = - 0 
xx yy 


(4.19) 


With this substitution equations (4.15) and (4.16) become 
identical. This happens because, to this order, the 
overall energy equation is satisfied trivially by the 
specification of isothermal and isochoric flow. 


51 



The solution for the translational temperature 


tensor becomes 





(td' 

Q 

H 

<M 



U J 

1 3 

3 

(3P' 


xD 

[ td' 

2pc 

V 

[dTj 

V 

3 

l 3 J 



0 

0 


(4.20) 


The deviation of the translational temperature from 
isotropy is determined by the dimensionless product xD, 
which expresses the ratio of the time required for 
thermal equilibration to the inverse of the strain rate. 
xD plays the same role in the energy relaxation equation 
as kD does in the molecular strain relaxation equation. 

The stress field required to maintain this flow can 
be obtained from equation (3.27). With the present 
approximations, in terms of the individual components we 
have for the shear stress 



and for the normal stresses 

a m - P e + Y P, - 3n 
xx k xx (A k 


- 3n 



£ 

xy 
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(4.21) 
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(4.22) 
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k yy [X k pj 


e + T P. O 
yy X k xy 


f- P, - 3 n P , 
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z z 


V ) £ 
xy xy 

(4.23) 

(4.24) 


Substituting the solutions obtained for 0 and £^, 


the shear stress becomes 

P k /"gPl ^ . A-B ( d „ „ „ 

a = - xD + — 7 — — P, - 3n P 

xy 2pc 3T 6 l X k p 

K v ^ ' v v 




(4.25) 


and the normal stresses become 
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and 


(A-B) 2 [" d 


z z 


12 


IT P u - 3n P 
14... k P 


(4.28) 


Since the shear rate is low, we have 


and 


ID << 1 
KD << 1 
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Thus the normal stresses can be seen to be of higher 
order than the shear stresses. 

For a dilute gas, because of the wide separation 
between molecules, the attractive forces become negli- 
gible, and 



The stress field then becomes simply 


a 

xy 


P 


TD 

3 


a 

XX 



a 


yy 



(4.29) 

(4.30) 

(4.31) 


a =0 (4.32) 

z z 

This is seen to be identical to the results obtained for 
a Maxwell fluid (Ref. 46). Though the strain rates 
required to generate them would be extremely high, these 
results indicate that normal stress effects can be 
expected in gases as well as in liquids. 

Normal stress effects in gases have also been cal- 
culated by Morgan and Kern (Ref. 14), who obtained off- 
diagonal terms of the stress tensor in axes coincident 
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with the principal axes of the strain rate tensor, and by 
Truesdell (Ref. 34), who did a kinetic theory analysis of 
Maxwell molecules. Reiner (Refs. 35, 36) claims to have 
measured normal stress effects in air, however, there is 
still some question concerning the interpretation of his 
results (Ref, 20), 

The result for the shear stress, equation C4.25), 
indicates that for slow shear rates the behavior of the 
fluid is Newtonian, with a viscosity 


V 


a 

= = 
D 


_i (jyp 

2pc 19T 
v '• 


P 



+ 


A-B [~d 

6 X 


3n 



(4.33) 


The first term represents the kinetic contribution and is 
proportional to the square root of the temperature. This 
is identical to the kinetic theory result for a gas com- 
posed of rigid elastic spheres. 

The second term is a combination of kinetic and 
intermolecular force contributions and decreases expo- 
nentially with temperature. For most liquids it is 
expected that the gas type contribution to the viscosity 
is small and that therefore the second term dominates the 
^ first . The viscosity can then be represented as 
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C4.34) 


an expression which is very similar in form to the hybrid 
equation for viscosity advanced by Macedo and Litovitz 
(Ref. 37). Their equation, 

y - A o 

where y is a correction factor required to correct for 
the overlap of free volume, was proposed to combine the 
desirable aspects of the free volume theory of Cohen and 
Turnbull (Ref. 38) with the reaction rate theory of 
Eyring (Ref. 39), both of which have proven to be very 
successful in predicting the behavior of different liquids. 
Considerable success has been achieved by the application 
of the Macedo-Litovit z equation to the temperature and 
pressure dependence of a wide variety of liquids, both 
organic and inorganic. The main difference between the 
Macedo-Litovitz equation and equation (4.34) arises in 
the pre-exponential factor, which for the Macedo-Litovitz 
equation is an adjustable parameter. A table, listing 
the values of and V Q / V £ required to fit the viscosity 
of many liquids, is presented by Macedo and Litovitz. 


ex 'e tf + Y ^ 
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From these values it can be seen that the exponential 
term is in general very large, and that for most liquids, 
therefore, k >> t. This justifies the neglect of the gas 
like term in the viscosity of liquids. 

It is important to note that the expression obtained 
for viscosity, particularly the pre-exponential term, is 
subject to all the previously enumerated limitations of 
cell theory and the equation of state. Further, it is 
necessary to account for the fact that the hard sphere 
volume is not a constant, but decreases with increasing 
temperature, since more energetic molecules can approach 
each other more closely than those with low energy. 

Equation (4.34) has been applied to two oils, for 
which the variation of viscosity and density have care- 
fully been measured over a wide range of pressure and 
temperature^ These are oils 17-D and 29-F of the ASME 
Pressure Viscosity Report (1953). The main reason for 
choosing these particular oils was because of their sim- 
ilarity to the oils used by Crook (Ref. 40) and Smith 
(Ref. 41) in their e las t ohydro dyn ami c lubrication exper- 
iments which are discussed in the next section. The mol- 
ecules of the oils are certainly not spherical, however, 
the presence of large proportions of cyclic hydrocarbons, 
such as the napthenes and the aromatics, could make this 
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a reasonable first approximation. J2?he success of the 
Cohen and Turnbull equation, which is a special case of 
equation (4.35) in describing the viscosity data of high 
molecular weight hydrocarbons to within the accuracy of 
the experimental data is certainly encouraging in this 
respect . 

The molecular weight of oil 29-F was not avail- 
able, consequently, it was estimated by the method of 
Hirschler (Ref. 42), to be approximately 440. The index 
n, for the variation of the potential energy with volume 
was taken equal to one, since this has been found to be 
a good approximation for many liquids (Ref. 3). It was 
further assumed that the cross -s t rains produced by the 
strain rate are small, and therefore A was put equal to 
one and B to zero . 

The energy barrier, E , and the close packed 

volume, v , were not known, and were therefore treated as 
o 

adjustable parameters. It has been found (Ref. 29) that 
E ^ is an increasing function of density, and that v q 
decreases with increasing temperature, accordingly, 
simple linear relationships displaying these characteris- 
tics were used. 

Gubbins and Tham (Ref. 29) have shown a reduced 
plot of the "activation energy", E^, as a function of the 
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volume for several liquids, and their curve, together 
with the points obtained for the two oils, is presented 
in figure (2). For the oils, the pour point parameters 
were used for the reduction of the data since the melting 
point parameters used by Gubbins and Tham were unavail- 
able. Nevertheless, the agreement of the values for the 
oils with those obtained for other liquids is quite good. 

The experimental and calculated variation of vis- 
cosity with temperature and pressure is shown for the two 
oils in figures (3) and (_4). The agreement is very good, 
except for high temperature and low pressures, which is 
precisely the region where cell theory cannot be expected 
to apply. 

4 . 2 Steady Shearing Flow at Arbitrary Strain Rates 


The analysis presented in the previous section is 
applicable to infinitesimally small strain rates. When 
the restriction, KD << 1, is removed, the equations are 
no longer linear. In this case we can expect that the 
stress will no longer be linear in the strain rate, and 
thus the flow will be non-Newtonian. 

Non-Newtonian behavior has been observed in liquid 
crystals (Ref. 43), colloidal suspensions (Ref. 44), 
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polymer melts and solutions (Ref. 45), and mineral oils 
in elas t ohy dro dynami c (.EHD) lubrication (Ref. 46). The 
theory, as developed so far, cannot be expected to hold 
for the first three cases, since for liquid crystals and 
polymers the molecules are extremely complicated and 
orientation effects play a large role, and for colloidal 
suspensions the flow of the particles making up the 
suspension must be considered (Ref. 47). There is hope, 
however, that the high shear rate flow of mineral oils 
can be described adequately by the present theory. This 
hope is reinforced by the success achieved in describing 
the oil viscosities over a wide range of temperature and 
pressure . 

The elas tohy drodynamic problem consists of the 
combination of the elastic deformation of the lubricating 
surfaces and the hydrodynamic behavior of the lubricant 
film. The shape of the film is determined by the defor- 
mation of the contacting surfaces, and since this film is 
generally very thin, the strain rates are extremely high. 
Figure (5) shows the shape and thickness of the lubricant 
film between two discs, measured by capacitance tech- 
niques (Ref. 40). The thickness is seen to be constant 
to within 10% in the region where the discs are deformed 
or in the Hertzian band. The ratio of the length of the 
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contact area to the film thickness is of order 100, and 
the width of the area is often 10 3 to IQ 4 times greater 
than the film thickness. The extremely high pressures 
encountered raise the viscosity of the oil in the contact 
zone by factors of 10 4 or more, thus reducing the escape 
of lubricant from the sides and front. These considera- 
tions indicate that the EHD flow of lubricants can be 
approximated quite well by plane Couette flow. 

To analyze this flow, we make the same assumptions 
as in the previous sections, except that we remove the 
restriction kD << 1. Also, since it was shown that the 
viscosity of the liquid is primarily determined by the 
molecular strain, the energy relaxation equations are not 
considered . 

The molecular strain relaxation equations now 

b e come 


for i = x , j = x 


Cl + e ) Cl + £ ) D e + e (1 + e ) — (e -e ) 

yy zz xy xy zz 2 yy xx 


(e - e ) + (e - e ) 

xx yy xx z z 


= 0 


( 4 . 36 ) 
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for i = y, j = y 
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xx zz xy xy zz 2 yy xx 
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for i = z , j = z 
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XX 
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(4.38) 


and for i = x, j = y 


3£ 

(1 + £ ) ( 2 + £ + £ ) 7 (£ - £ ) + ^ D . 

zz xx yy 4 xx yy < 2 

(4.39) 


Substituting the expression obtained for from 

equation (4.38) into equation (4.36) and (4.37) makes 
the latter equations identical for the same reason which 
was discussed in the previous section, therefore we again 
use equation (2.3). 

With some algebraic manipulations, the equations to 
be solved can be written as 
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(A. 42) 


(4.43) 


In these expressions, K is still a function of the mole- 
cular strain, and this relationship is yet to be deter- 
mined . 


Gubbins and Tham (Ref. 29) give the following 
description of the jump of a molecule into a neighboring 
hole. In figure (6), for the molecule at 1 to jump to 3, 
it is necessary to overcome the energy barrier at 2 
caused by the steep repulsive potentials of the neighbor- 
ing molecules. This barrier is clearly a function of the 
molecular spacing in the directions normal to the direc- 
tion of the jump. Thus for a jump to occur in the z 
direction, the energy required should be a function of 
tne spacing in the x and y directions. A plausible 
assumption would be that depends on the area of the 
gap as seen by the jumping molecule and thus the varia- 


tion of can be expressed as a function of the product 


A A . In principal axes, a, we can then write for small 
x y 


strains that 
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(4.44) 
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^ e v 

a a " E ve“ 2 d& £ aa 
e 

where dE / dSL can be determined from figure (2) . 
v e 

For pure shear, a sphere is distorted into an 
ellipsoid whose axes correspond to the principal axes of 
strain. We therefore assume that the energy required for 
a jump in a particular direction is also distributed 
ellipsoidally . Hence we can write 


E y (.<}>, <10 


E ve + 


i dE 
e v 

2 dSL 
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I, 




sin 2 <j>cos 2 \p sin 2 <J>sin 2 i p 

- - + - - + 

[ Cl + Ei) 2 (1 + e 2 ) 2 



(4.45) 


* * * 

where Ei, e 2 , and £3 are the principal strains obtained 
from the roots of the equation 


£ 3 £ 2 £ 

e - I £ e - II e e - IIl e = 0 . (4.46) 


The invariants of the tensor (.E^ + S^j) are 
I e = E ii + 5 ii 

' I !(e li + 6 ii )U 33 + S Ja >- Ce ij + + 

and 


III £ = det(e. + 6 ) . 
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This expression for £(<}>,$) can be used in equation (3.48) 
to integrate for the jump probability over all direc- 
tions. Similar arguments can be used to determine the 
change in the probability of finding sufficient free 
volume for a jump, as a function of direction, and the 
same ellipsoidal distribution is used for v q / v f (<p , ip ) in 
equation (3.48) . 

Equations (4.40-43) can be solved numerically to 

obtain the strains as a function of < D, where K is the 

e e 

molecular strain relaxation time for the unstrained 
liquid . 

The solutions are obtained by the following proce- 

_ 4 

dure. First, a small value for ic^D, typically 10 , is 
chosen, and the strains are computed with the low strain 
rate assumption, equation (4.12). These strains are used 
as a first approximation in solving equations (4.40-43) 
by a Newton-Raphson technique for systems, where KD is 
assumed constant. With the new strains, a new KD is 
calculated from equations (3.46) and (3.48) and compared 
with the previous kD. If the two values for KD are close 
it is assumed that the method has converged to the 
required solution, if not, the new value for KD is used 
to obtain another approximation to the strains. This 
procedure can be continued until convergence is obtained. 
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When a higher K q D is chosen next, the previous solution 
is used as the first guess for the strains. 

Once the required strains have been obtained, the 
stress field can be calculated from the following 
exp re s s ions . 
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(4.48) 

(4.49) 

(4.50) 


It is convenient to non-dimensionalize the stresses by 
dividing them by the expression 



Arifi P 

6 [X k 



(4.51) 


where is later found to be the limiting high frequency 
shear rigidity. For a Newtonian fluid the dimensionless 
shear stress is just given by 



(4.52) 
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The numerical solutions obtained for the variation 
of shear stress with shear rate are compared with 
Dyson's (Ref. 46) interpretation of the EHD lubrication 
experiments performed by Smith (Ref. 41) and Crook 
1 (Ref. 40) . 

Smith's experimental apparatus consisted of a motor 
driven cylindrical roller in contact with a spherical 
roller which was allowed to rotate freely in ball 
bearings. The axes of the two rollers could be skewed 
relative to each other, thus varying the sliding speed 
between them. Measurements of the axial force on the 
spherical roller were used to determine the frictional 
force. Heating elements inside the housings of the 
rollers, and a thermocouple in the air space near the 
point of contact of the rollers were used to maintain and 
regulate the temperature of the system. The film thick- 
nesses were not measured, however, Dyson calculated them 
using the approximate theory of Grubin (Ref. 48). He 
estimates that this theory should be sufficiently 
accurate to give the film thicknesses to within a' factor 
o f three . 

Dyson further introduced corrections to the exper- 
imental points for temperature variations in the contact 
zone, determining what the results would have been if the 
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conditions had been isothermal. 

Crook performed his experiments by loading together 
pairs of discs, and measuring the frictional force as 
the sliding speed was varied. The film thicknesses were 
determined from the measured capacitance between the 
discs, and the surface temperatures were measured with 
thermocouples. Dyson also corrected Crook’s results for 
thermal effects. 

There is yet another correction which must be made 
before the calculated and experimental results can be 
compared, and that is to account for the variation of 
pressure in the contact zone. 

It can be shown that under high loads, the pressure 
distribution in the contact region during EHD lubrication 
is very close to that which would exist for static 
loading with no lubricating film. This "Hertzian" pres- 
sure distribution is given by 


P 




( 4 . 53 ) 


where b is the half width of the contact region. Using 
Grubin's assumption that the viscosity variation with 
pressure is nearly exponential, 
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(4.54) 


y = y exp(BP) 


where B depends only on temperature, it is possible to 
develop a relationship between the mean effective vis- 
cosity, and the viscosity at the maximum Hertzian pres- 
sure. The ratio between the two can be expressed as 


( : 


exp <3? 


fi -^1 




mean 


max 



(4.55) 


exp (BP ) 
m 


The experimental results were all adjusted to the condi- 
tions prevailing at the maximum Hertzian pressure, which 
were also the conditions for which the calculations were 
performed . 

The properties of the mineral oils studied by Smith 
and Crook were not available, therefore the oils showing 
the greatest s imi lar i ty to the actual ones used were 
chosen from the ASMS Pressure Viscosity Report (1953) . 

Oil 17-D of the report was chosen to be representative of 
Smith's oil, while Crock states that oil 29-F was most 
s irnilar to his . 

The experimental conditions were all out of the 
range of the ASME data, however, some cases were close 
enough for a reasonable extrapolation to be attempted. 

For this reason some of Smith's high pressure data could 
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not be compared with the theory. The most difficult 
quantity to find an extrapolated value for is 

Z dE 

^ = 2kT dZ (A. 56) 

\ • • 

which is required to evaluate the change in the jump 

probability due to the molecular strain. From figure (2) 

we see that E rises very rapidly as V/V is decreased, 
v m 

consequently it is extremely difficult to extrapolate to 

small values of V/V . This problem is compounded further 

m 

by attempting to determine an extrapolated value for the 
slope of the curve. For this reason, f is treated as an 
adjustable parameter in the calculations. The values for 
the hard sphere diameters and close packed volumes are 
obtained from the relationships used to fit the low shear 

rate viscosity data. The conditions of the experiments, 

. / 

and the quantities used for the calculations are shown in 
Table I. 

The experimental results obtained by Smith and Crook, 

for the variation of shear stress as a function of strain 

rate are shown plotted in dimensionless form in figure 

(7). Since the theory shows that in addition to being a 

function of the dimensionless strain rate, K D, the 

e 

stress also depends upon the dimensionless quantities 
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it is surprising that the 


¥, d/X, E /kT. and P. /3nP , 
v k p 

results appear, within the scatter of the data, to fall 
on a single curve. This is particularly unexpected in 
view of the dissimilarity of the two oils used in the 
experiments and the fact that yet two other oils were 
used for the reduction) of the data. This result may be 
fortuitous since the data plotted in dimensional form is 
originally close to a single curve, however, reduction 
to dimensionless form reduces the range of scatter from 
a factor of three to a factor of two, or roughly thirty 
per cent. 

The general features of Smith's results may be seen 
in figure (8) where the data obtained using steel rollers 
is plotted in dimensional form. An increase in the pres- 
sure at constant temperature or a decrease in the temper- 
ature at constant pressure raises the stress level at 
which departure from Newtonian flow occurs. The behavior 
of the stress level is similar to that of G , therefore 
non-dimens ionaliz ing the stresses with can be expected 
to bring the results together, though there is no a priori 
reason for expecting them to collapse onto a single curve. 

The experimental conditions for all but the lowest 
pressure and highest temperature were too far out of the 
range of the ASME data to permit extrapolation with any 
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degree of confidence, accordingly, only this case is 
shown in dimensionless form in figure (9). The calcu- 
lated curve, together with the value of ¥ required to 
match the data is also displayed, and the fit is seen 
to be quite good. 

In figure (10) Smith's results for a spherical alu- 
minum free roller rolling on a cylindrical steel driven 
roller are shown. These results indicate that the shear 
stress at which departure from Newtonian flow occurs is 
independent of the rolling speed, even though, as has 
been discussed by Dyson, the rolling speed does affect 
the apparent viscosity. The reason for this behavior is 
still not understood. The calculated curve can again be 
fitted quite well to the measured data. 

Crook's results, together with calculated curves for 
his average conditions, are shown in figure (11). Since 
these measurements were taken over a much narrower range 
of strain rates than Smith's, the effect of the scatter 
is much more pronounced. 

In figure (12) we show the calculated results for the 
change in shear stress when some of the non-dimensional 
parameters are varied. A decrease in the values of ¥, 
E^/kT and d/X produce an increase in the dimensionless 
stress at which departure from Newtonian flow takes place, 
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however, the variation of E^/kT and d/X, throughout the 
entire range encountered in the experiments, is insuffi- 
cient to shift the calculated results out of the range of 
scatter of the data. This indicates that, in this range, 
the stresses are fairly insensitive to changes in E^/kT 
and d/X, and that the largest variation occurs as a 
result of changes in 't'. Altering the ratio of the kinetic 
to potential pressure is found to have a negligible effect 
on the calculated results. The value of Y required to 
match the experimental data is a strong function of the 
values of G used to non-dimens ionali ze the stresses, 
thus any inaccuracies in G m produce a corresponding error 
in ¥. 


Neither Crook nor Smith attempted to measure normal 

stress effects, nevertheless, since the theory predicts 

their presence, the calculated results are indicated in 

figure (13). In general, we see that a and a are 

yy xx 

both negative, suggesting that the flow causes an increase 

in the mean compressive stress. Further, since 

(o - a ) > 0, the fluid is in tension along stream- 
xx yy ° 

lines and in compression normal to the stream lines along 


the velocity gradient, when compared to the mean normal 

stress. Since a is a factor of approximately a hundred 

z z 

smaller than either a or C , the relative magnitude of 

xx yy 
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the stresses are very similar to those obtained for a 
Maxwell fluid. This resemblance to a Maxwell fluid is 
also shown by the low shear rate results, where it is 
seen that the normal stresses depend on the square of the 
strain rate. Changes in E^/kT have very little effect on 
the normal stresses, however, decreases in V and d/X lead 
to an increase in their magnitude at the corresponding 
strain rates. 

The variation of molecular strain with strain rate 

is shown in figure (14). The magnitude of the strains, 

-2 

for all the cases considered, is always less than 10 

indicating that the small strain assumption should be 

valid for an entire range of flow. The stress, however, 

is a function of V.. = 6 .. + (1 + d/X)e , which for 

ij ij ij 

large values of d/X can become large. Thus it is neces- 
sary to evaluate e without making the small strain 
assumption . 

We can conclude that the present theory is capable 
of describing non-Newtonian behavior in a hard sphere 
fluid. In this theory the decrease in apparent viscosity 
with increasing strain rate can be caused by two effects. 
One is the increase in the jump probability of a molecule 
in a strained fluid, which takes place at high densities 
and high strain rates. This mechanism is characterized 
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by a high, value of 

The other is that as the strain rate increases, the 
strain eventually reaches a limiting value. When this 
value is reached, an increase in the strain rate no 
longer causes an increase in the potential term of the 
stress tensor. This mechanism, as can be seen from 
equation (4.47) is characterized by a low value of 
(d/X) (P./3nP ) . 

k p 

There is yet another mechanism which could be 
responsible for non-Newtonian behavior of liquids, but 
since it entails the addition of some more assumptions, 
it has not been included in the present theory. In 
associated liquids it is possible that the rate of making 
and breaking of hydrogen bonds is a function of the mol- 
ecular strain or the strain rate. Since the jump prob- 
ability of a molecule is dependent on whether or not it 
is hydrogen bonded to others, a decrease in the number of 
bonded molecules can lead to a decrease in the apparent 
viscosity of the liquid. Similar arguments could very 
well apply to any liquid whose molecules form clusters. 
The sizes and the number of clusters could be dependent 
on the strain rate,* and -since a molecule in a cluster 
would probably have a smaller jump probability than a 
free molecule, a decrease in the size and number of 
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clusters would lead to a decrease in the apparent 
viscosity . 


4 . 3 Oscillatory Shear Flow 


Consider the propagation of plane oscillating shear 
waves in a liquid. Let x be the direction of motion of 
the liquid and y be the direction of propagation of the 
shear wave. Assume that the amplitudes of the oscilla- 
tions are sufficiently small that products of perturba- 
tion quantities may be neglected. Then the equations 
describing the flow can be linearized to give: 


the molecular strain relaxation equation 

3 £ 3£ . _ 3u 

xy xy A-B x 

3t K ~ 2 3y 

the energy relaxation equation for no conduction 


(4.5 7) 
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(4.58) 


the momentum equation 
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(4.59) 


and the stress 
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A K p xy 


(4.60) 
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Assume that any perturbation quantity Q, varies as 


Q = Q exp(io)t + m*y) 


(4.61) 


where u is the angular frequency and m*, is defined as 

iw 


m* = - a - 


(4.62) 


Here a is the absorption coefficient and c is the 
propagation velocity. It is further convenient to let 


and 


k’ = K/3 
x ' = T/3 


(4.63) 

(4.64) 


The equations which must be solved then become 


iue + e /<' = m*u 

xy xy 2 x 
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(4.68) 


Solving equations (4.65) and (4.66), we find that 
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; _ A-B X 

'xy 2 It iuK' 


(4.69) 
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and 


6 

xy 

Substituting the expression for the stress, (4.68), in 
the momentum equation (4.67), and then using the expres- 
sion obtained above for e and 9 we get 

xy xy 
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(4.71) 


In practice the quantity measured in oscillating shear 


experiments is the shear mechanical impedance, Z g , 
defined by 


Z = — = R + iX . (4.72) 

s u s s 

X 


From the momentum equation we get 



9 


(4.73) 


which can be evaluated in terms of the fluid properties 
by using equation (4.71). Another commonly used quantity 
is the shear modulus, G, which is obtained from the 
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equation of motion for an elastic solid, 


P 


if! 

8t 2 



(4. 74) 


where £ is the displacement in the x direction. To 
account for dissipation, G is made complex, and solving 
for oscillatory motion, equation (.4.74) gives the complex 
shear modulus 



(4.75) 


which can also be obtained from equation (4.71). The 
real and imaginary parts of the shear mechanical imped- 
ance can then be expressed as 


and 



(4. 76) 


(4.77) 


For a Newtonian fluid, the equation of motion is 


9u 8 2 u 

p nr = M nr 2 

and for oscillatory motion this yields 


(4.78) 
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(4.79) 


( iti)| 2 iuy 

l a *j P 


Comparing this expression with equation (4.75) we see 
that for a Newtonian fluid 

G* = iuy (4.80) 

which is purely imaginary. 

For a loss-free solid. 


G* = pc' 


(4.81) 


which is real. Here c g is the propagation velocity of 
the shear wave. 

According to the present theory, the complex shea 
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Taking the low frequency limit, we find that 


lim G* 
to -* 0 


= iw 


3 _!k 
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and comparing this to equation (4.33) we see that 


lim G* 
00 -*• 0 


iooy 


(4.85) 


This expression is the same as equation (4.80), thus in 
the low frequency limit our theory predicts Newtonian 
behavior. 

In the limit of high frequency, we get 


lim G* 
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which is real, and independent of frequency. Thus in 
this limit, our theory predicts that the fluid behaves 
like an elastic solid. 

The highest frequencies which can presently be 

obtained by experiment are in the neighborhood of I GHz, 

_ 1 3 

and since T ' is oforder 10 seconds, the terms 
involving the product cot' are very small. The molecular 
strain relaxation times of many liquids are of order 10 8 
or greater, thus these relaxation effects can be measured 
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The measurable complex shear modulus is 


G * = P k ' 3nP 

2 A k 
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which describes exactly the behavior of a single relaxa- 
tion time Maxwell fluid whose infinite frequency shear 
modulus is 


G 

OO 


A-B 

2 



- 3nP 

Pj 


(4.88) 


With this result, the low frequency, or low strain rate, 
viscosity can be expressed as 




G t' 


OO 


(4. 89) 


The very high frequencies required to study molecular 
relaxation processes can best be obtained by utilizing 
piezoelectric crystal transducers. Mason et at. (Ref. 49) 
developed a technique for obtaining the shear mechanical 
impedance of a liquid by measuring the reflection coeffi- 
cient of an incident shear wave at the crystal-liquid 
interface. This method has been employed and modified by 
many others to perform experiments on a wide variety of 
liquids (Ref. 9). 

So far only two liquids have been found which conform 
to the single relaxation time behavior predicted by both 
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the present theory and the Maxwell fluid model. These 
are molten zinc chloride (Ref. 50) and boric oxide 
(Ref. 51). The measurements on boric oxide were carried 
out over a wide range of temperatures, and it was found 
that above 800°C the fluid exhibits single relaxation 
time behavior, while below this temperature the width of 
the relaxation region increases with decreasing tempera- 
ture. 

Up to the present time, all other liquids, on which 
oscillatory shear experiments have been performed, 
exhibit a relaxation region which is broader than that 
predicted for a single relaxation time process. This 
behavior has generally been analyzed by postulating that 
the viscoelastic behavior can be represented by the summa 
tion of discrete single relaxation time processes, and 
then determining the spectrum of relaxation times 
required to fit the data. As has been pointed out by 
Barlow et al . (Ref. 52), this procedure does not yield a 
unique distribution, nor does it give a physical descrip- 
tion of the process. As an alternative approach, Barlow, 
Erginsav and Lamb (Ref. 53) have suggested that the shear 
me chan i-c a 1 --impedance of a visco e lastic liquid can b e -- r . - 
represented by a parallel combination of impedances 
characterizing a solid and a Newtonian fluid. This 
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procedure has had a fair degree of success in describing 
the relaxation behavior of some fluids. A comparison of 
the results obtained, for the real and imaginary parts of 
the shear mechanical impedance, using this model and the 
single relaxation time theory is presented in figure (15) 

An attempt at a physical explanation for the width 
of the relaxation region has been provided by Litovitz 
and McDuffie (Ref. 54) who in a qualitative way attribute 
it to the requirement of cooperative motion of neighbor- 
ing molecules to provide the free volume necessary for 
molecular rearrangement. 

Another possible explanation for this behavior could 
be the effect of molecular clusters in the liquid. At 
high temperatures the number and sizes of the clusters 
would be small, meaning that all the molecules would 
essentially be free. This would lead to a single relax- 
ation time behavior. As the temperature is decreased, 
the number and sizes of clusters would increase, making 
the liquid a mixture consisting of a wide range of 
cluster sizes. This would lead to a spread in the range 
of the relaxation region with decreasing temperature, as 
observed in B2O3. 

There does not yet appear to be a satisfactory quan- 
titative physical explanation of the width of the visco- 
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elastic relaxation region. 

Barlow and Lamb (Ref. 55) have experimentally 
studied the behavior of three lubricating oils possessing 
high, medium and low viscosity indices. All of these 
exhibited a relaxation region which was broader than that 
predicted by single relaxation time theory, however, it 
is of interest to see how well the present theory pre- 
dicts the magnitude of the limiting high frequency shear 
modulus, G^. The viscosities, at 30°C, of the high, 
medium and low viscosity index oils were respectively 
3.54, 2.08 and 4.52 poise, while those of oils 29-F and 
17-D were 1.6 and 0.12 poise. Thus we would expect oil 
29-F to be more similar to the oils used by Barlow and 
Lamb than oil 17-D, and thus provide a better comparison. 

Using the values of d/A obtained by fitting the low 
shear rate viscosity data, the calculated results for 
were 7.2 x 10 9 and 3.9 x 10 9 dyne/cm 2 for oils 29-F and 
17-D respectively. The measured values obtained by 
Barlow and Lamb were 7.0 x 10 9 , 7.8 x 10 9 and 7.0 x 10 9 
dyne/cm 2 for the high, medium and low viscosity index 
oils in that order. Since there were no adjustable 
parameters in the calculation of G , these results must 
certainly be considered to be encouraging. 

Comparing the experimental results for oscillatory 
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and steady shear flows it has. been noted (Ref. 46) that 
departure from Newtonian flow occurs at much lower strain 
rates for the latter case than the former. This feature 
cannot be explained on the basis of the Maxwell model 
since the values of required to describe the two flows 
by this model differ by a factor of about 10 3 . The pres- 
ent theory, however, can account for this disparate 
behavior on the basis of the difference in magnitude of 
the molecular strains encountered. One consequence of 
this is that the behavior of a liquid under conditions of 
steady shearing flow cannot, on the basis of the present 
theory, be deduced from oscillatory shear experiments. 

4 . 4 Absorption and Dispersion of Longitudinal Waves 

Consider the propagation, in the x-direction, of 
plane longitudinal periodic waves of infinitesimal ampli- 
tude. In most liquids, excluding the liquid metals whose 
thermal conductivities are high, the effects of conduc- 
tion are very small compared to the effects of viscosity 
(Ref. 5). Further, when the period of the sound wave is 
large compared to the relaxation times of the fluid, the 
effects of viscosity and thermal conductivity on the 
absorption add linearly. Therefore, to simplify the 
analysis, thermal conduction is ignored. 
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The equations which must be solved are then: 
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momentum 
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and the dynamical equation of state 
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We now write the local equilibrium quantities as 


T = T (1 + 0) (4.96) 

o 

and 

p = p Q Cl + s ) (4.97) 

where the subscript o refers to the mean values, and 0 

and s are much less than one. The derivative of S with 

xx 

respect to x, can then be written in terms of the deriv- 
atives of the perturbation quantities, and the momentum 
equation can be expressed as 
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If we now let the perturbation quantities vary as 
Q = Q exp (m*x + iut) 


and make the substitutions < ' = k/ 3 and T ? = t/3, the 
algebraic equations which must be solved are: 
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The solutions for the perturbation quantities are 
found to be 
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Substituting these expressions into equation 
(4.103), and separating the real and imaginary parts, the 
complex longitudinal modulus 
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can be expressed in terms of physical quantities as 
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The low frequency sound speed ts 
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Using standard thermodynamic relations this can be 
rewritten in the familiar form 



(4.112) 


where y is the ratio of specific heats. 
If we make the approximation 
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a condition fulfilled by most liquids (Ref. 9), the 
propagation velocity can be expressed as 


c 2 = 


M* 


(4.113) 


and looking at equation . (4 . 109) we see that velocity 
dispersion is predicted. Since at the present time 
measurements cannot be performed at frequencies of com- 


parable magnitude to the collision frequencies in liquids, 
the terms involving wt are small and the measurable 
velocity dispersion is given by 
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which using equation C4.88) can be rewritten as 
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This is the result for a viscoelastic fluid with no bulk 
viscosity (Ref. 56). 

The absorption, a, can be obtained from the relation 
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and it can be seen to exhibit two maxima, one at wic’ = 1 
and the other at WT ’ *» 1. For low frequencies the absorp- 
tion becomes 


a 



(4.117) 


and comparing this to equation (4.35) we find that 
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P c 3 
o o 


(4.118) 


which is just the classical result for the case of no 
heat conduction. We can thus see that the present theory 
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predicts a zero bulk, or volume viscosity for a hard 
sphere liquid. This disagrees with experimental results 
which indicate that monatomic liquids such as liquid 
argon (Ref. 57) and several liquid metals (Refs. 58, 59) 
exhibit excess absorption. The reason for this discrep- 
ancy is that all the relaxation mechanisms have not been 
considered in the theory so far. 

Generally in nonassociated liquids the major causes 
of excess absorption are rotational isomerism and vibra- 
tional relaxation, whereas in associated liquids the 
effects of structural relaxation are most important. 
Structural relaxation has also been suggested to be 
responsible for the existence of volume viscosity in 
monatomic liquids (Refs. 9, 58, 59). This effect could 
arise in the following fashion. 

It is well known, through the property of critical 
opalescence, that liquids near the critical point exhibit 
large, local density fluctuations caused by the continuous 
agglomeration and breakup of large clusters of molecules. 
These clusters often contain on the order of millions of 
molecules. Away from the critical point, clustering 
should still occur to some extent in most liquids. The 
equilibrium numbers and sizes of the clusters are deter- 
mined by the local thermodynamic conditions, and any 
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changes in these conditions lead to potential energy 
changes due to the varying structure of the liquid. There 
is therefore associated with this process a characteristic 
time, or structural relaxation time, which describes the 
equilibration of the cluster size and number distribution 
after a perturbation of the equilibrium conditions. This 
process can therefore lead to absorption in excess of 
that predicted by the classical mechanisms of viscosity 
and conduction. 
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CHAPTER V 


SUMMARY AND CONCLUSIONS 

A relatively simple theory to account for the dynam- 
ical behavior of liquids has been proposed, using a mole- 
cular viewpoint, to obtain general dynamical equations 
similar to the phenomenological relaxation equations 
first advanced by Maxwell (Ref. 10). The theory is based 
on the concept that in a continuously straining liquid 
distortions from spherical symmetry of the average mole- 
cular spacing, and the average random translational 
energy of the molecules, take place. This leads to the 
definition of a molecular strain tensor to describe the 
deviation of the radial distribution function from its 
equilibrium value, and a translational temperature tensor 
to describe the directional dependence of the random 
translational energy of the molecules. If all disturb- 
ances are removed from the fluid, the radial and velocity 
distribution functions eventually revert to spherical 
symmetry at rates characterized by the relaxation times 
for each process. The relaxation time for equilibration 
of translational kinetic energy has been identified with 
the collision time, while the time required to obtain 
isotropicity of molecular spacing has been related to the 
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time required for a molecule to jump into a neighboring 
vacancy . 

A dynamical equation of state has been obtained, 
using statistical mechanical techniques, by postulating 
that directionally dependent, quasi equilibrium, velocity 
and radial distribution functions could be written. The 
stress tensor has thus been related to the molecular 
strain and translational temperature tensors by evalu- 
ating the partition function obtained from simple cell 
theory and considering the distortion of the free volume. 

The equations obtained from the proposed theory have 
been solved for the problems of steady shearing flow, the 
propagation of oscillating shear waves, and the propaga- 
tion of longitudinal waves. 

For steady shear flow the solutions indicate that at 
slow shear rates the flow is Newtonian, and an expression 
for the viscosity of a liquid similar to the very success- 
ful Macedo-Litovit z equation is ybtained. This equation 
has been applied, with good success, to the pressure and 
temperature dependence of the viscosity of two lubricating 
oils . 

At high shear rates, the decrease in apparent vis- 
cosity with increasing strain rate observed by both Smith 
and Crook (Refs. 40, 41) in their EHD lubrication experi- 
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ments, can be described using the present theory. The 
reduction in apparent viscosity can be attributed to the 
enhanced jump probability of the molecules caused by the 
distortion of the radial distribution function from 
spherical symmetry. Calculations also predict the exist- 
ence of normal stresses whose relative magnitudes are 
similar to those obtained for a simple Maxwell fluid. 

Under conditions of oscillatory shear, results 
obtained from the present theory indicate that a liquid 
behaves like a Newtonian fluid at low frequency and like 
an elastic solid at high frequency. The relaxation 
region between these limits is described by precisely the 
same expression as for a Maxwell fluid with a single 
relaxation time. Experiments indicate, however, that for 
most liquids the relaxation region is in general broader 
than that described by single relaxation time theory. 

This discrepancy has been attributed to the presence of 
molecular clusters. The calculated numerical values of 
the high frequency shear modulus for the two lubricating 
oils considered in the high shear rate experiments, turn 
out to be of the same order of magnitude as* those meas- 
ured by Barlow and Lamb (Ref. 55) in oscillatory shear. 

The theory which has been proposed also provides a 
satisfactory qualitative description of the propagation 
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of longitudinal waves in a liquid. At low frequencies 
both the sound speed and the absorption are found to have 
classical values, and at higher frequencies velocity dis- 
persion with a single relaxation time is obtained. The 
experimentally observed excess absorption in monatomic 
liquids, indicating a non-zero bulk viscosity, is again 
attributable to the presence of molecular clusters. 

A general feature of all the cases studied so far is 
that classical results are predicted for low strain rates, 
and that as the strain rate is increased, the behavior of 
the liquid becomes viscoelastic. The present theory is, 
however, superior to the phenomenological viscoelastic 
theories in that the elastic moduli are all expressed in 
terms of molecular and thermodynamic properties. Thus 
from a qualitative viewpoint a good description of the 
dynamical behavior of liquids is provided. Nevertheless, 
some discrepancies in the quantitative results exist. 

These are due to a large extent to the fact that other 
relaxation processes, such as thermal and structural 
relaxation, have not been included. 

Thermal relaxation effects can be incorporated into 
the present theory in precisely the same way as they have 

r 

been added to the classical theory as described by 
Herzfeld and Litovitz (Ref. 5). This is done by adding 
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the applicable energy mode to the energy equation and by 
writing relaxation equations connecting this mod': to the 
other modes between which energy exchange occurs. The 
results which have been obtained previously by relaxation 
theory are then also applicable to the present theory. 

Structural relaxation is at present not nearly as 
well understood as thermal relaxation. One promising 
approach in trying to gain some understanding of this 
process is to view the liquid as a mixture of molecular 
clusters of varying density, and defining a distribution 
function which describes the probability that a molecule 
is located in a cluster with a particular mean spacing. 
This method has been used with some success by Kerley 
(Ref. 60) to describe the thermodynamic properties of 
liquid argon and hydrogen. If it would be possible to 
relate this distribution function to flow conditions, it 
should be feasible to write relaxation equations which 
describe the behavior of the clusters under dynamic 
conditions, and in such a manner account for the excess 
absorption and the width of the relaxation region. 

Further improvements which -could be made to the 
theory are to use a better approximation to the partition 
function than that provided by simple cell theory and to 
obtain a more accurate expression for the collision 
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frequency in the liquid. Experimentally it would be 
useful if high shear rate and oscillating shear experi 
ments could be performed on the same liquid, and if 
thermodynamic data and transport properties for this 
liquid could be obtained in the high pressure range of 
the steady shear flow experiments. 
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Figure 2. The energy required for a molecular jump as a 
function of volume 
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VISCOSITY (poise) 



Figure 3. The viscosity of 0il-17D as a function of 
temperature and pressure 
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Figure 4. The viscosity of Oil 29-F as a function 
of temperature and pressure 
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Figure 5. Typical film thickness variation 
in EHD lubrication 



Figure 6. Mechanism of a molecular jump 
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Figure 7. Measured variation of the dimensionless shear stress as 



(dyn/cm 2 ) 

Figure 8. Variation of shear stress as a function of strain rate and load 
obtained by Smith using steel rollers 
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Figure 9. Comparison of measured and calculated shear stress variation as 
function of strain rate. Results of Smith using steel rollers. 
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function of strain rate. Results of Smith using steel on hard aluminum 
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Figure 11. Comparison of measured and calculated shear stress variation as 
function of strain rate. Results of Crook. 
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Figure 12. The effects of varying parameters on the calculated shear stress as 
function of strain rate 
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Figure 14 . 


Calculated variations of molecular strain 
in steady shear as a function of the 
strain rate 
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Figure 15. Real and imaginary components of the shear impedance in 
oscillatory shear 
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